Signatures of modulated pair interaction in cuprate superconductors 
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Abstract 

Recent low-temperature scanning tunnelling spectroscopy experiments on the surface of Bi2Sr2CaCu2 08+i: have 
revealed a strong positive correlation between the position of localized resonances at -960 meV identified with 
interstitial oxygen dopants and the size of the local spectral gap. We review efforts to understand these correlations 
within a model where the dopants modulate the pair interaction on an atomic scale. We provide further evidence 
for this model by comparing the correlations between the dopants and the local density of states with experimental 
results. 
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1. Introduction 

High-resolution scanning tunnelling (STS) spec- 
troscopy on BSCCO-2212 has provided some of the 
most provocative data on cuprate superconductors 
in recent years, including local information on impu- 
rity states, charge ordering, quasiparticles in the vor- 
tex state, and nanoscale inhomogeneity. In partic- 
ular the observation of nanoscale inhomogeneity in 
the electronic structure of BSCCO-2212 has raised 
great interest due to its possible relevance to under- 
standing how the cuprates evolve from the doped 
Mott insulator to the superconducting state. Re- 
cently, STS provided intriguing results which may 
bear on this question, by the observation of random 
localized resonances at a bias of -960 meV which 
were argued to be states associated with interstitial 
oxygen atoms[l]. These defects dope the BSCCO- 
2212 crystals, along with a rather complex mix of 
other dopants[2]. From these experiments a strong 
and positive correlation between the magnitude of 
the gap and the O dopant positions has been in- 



ferred by McElroy et al. [1]. This result was unantic- 
ipated, since it liad been suggested [3] that each O 
might overdope the crystal locally by providing two 
holes, corresponding to a smaller local gap in anal- 
ogy with the average gap decrease with overdoping 
known from tunnelling experiments [4]. Additional 
salient characteristics of the STS-spcctra arc their 
remarkable particle-hole symmetry and their homo- 
geneity at low energies. 

Nunner et al. [5] subsequently argued that these 
main characteristics of the STS-spectra cannot be 
explained by a model where the dopant atoms act 
simply as potential scatterers. Instead they sug- 
gested that the hypothesis of enhanced pairing 
around each 0-dopant can explain the observed 
correlations, in particular the negative correlation 
of the coherence peak height with the magnitude of 
the gap, i.e. the fact that the coherence peaks are 
suppressed in regions of large gap and enhanced in 
regions of small gap. The underlying physics within 
this model is Andreev-scattcring of the quasiparti- 
cles at order parameter modulations. In regions of 



Preprint submitted to Elsevier Science 



6 February 2008 



suppressed order parameter a new resonance forms 
below the gap edge [5] which has been used to ex- 
plain the extremely sharp coherence peaks observed 
by Fang et al. [6] in small gap regions. In regions 
of enhanced order parameter exactly the opposite 
happens and the coherence peaks are strongly sup- 
pressed [5] . The fact that the scattering mechanism 
is of Andreev type can also account for the remark- 
able particle-hole symmetry of the STS-spectra and 
for the observed homogeneity at low energies since 
nodal quasiparticles are less affected by order pa- 
rameter variations than antinodal quasiparticles. 
In addition, modulation of the charge is found to 
be quite weak, in agreement with experiment [1]. 
None of these effects can be found in models of 
conventional potential disorder. 

Further comparison of the pair disorder model 
with experiment revealed a natural explanation of 
the relative strengths of the quasiparticle interfer- 
ence peaks observed in Fourier transform scanning 
tunnelling spectroscopy (FT-STS) [7] and showed 
that the broadening of thermodynamic transitions 
induced by the gap modulations observed by STS 
was compatible with specific heat measurements on 
BSCCO-2212 [8]. This serves as an additional ev- 
idence that the nanoscale inhomogeneity observed 
by STS is indeed characteristic of the bulk in the 
BSCCO-2212 system. 

The success of this hypothesis raises the question 
of whether one can learn about the origins of pair- 
ing by studying what modulates it directly [9]. If 
one were to know the changes in electronic structure 
induced by a dopant or other impurity, and could 
establish a direct link to the change in the super- 
conducting order parameter, one would have valu- 
able information with which to constrain theories 
of superconductivity. A possible microscopic origin 
of the modulated pair interaction might be the dis- 
tortion of the crystal lattice around each O intersti- 
tial. These local lattice distortions might change lo- 
cal electron-phonon interaction matrix elements or 
modulate the nearest neighbor hopping t and as a 
consequence the superexchange constant J between 
nearest neighbor Cu-atoms (in a one-band model 
e.g. J — At'^/U). Within a strong-coupling picture 
one might assume that the pair interaction is pro- 
vided by the superexchange J and any modulations 
of J would directly translate into a modulated pair 
interaction. This scenario has been further worked 
out by Zhu [10] with results similar to the ones found 
by Nunner et al. [5] . 

To understand the detailed changes in electronic 



structure around a dopant atom, one can perform 
an ab initio density functional theory (DFT) calcu- 
lation of the structure with and without the inter- 
stitial present. A first attempt of this type for an 
oxygen interstitial in the BSCCO-2212 system was 
made by He et al. [11], who found that the lowest 
energy location for such a dopant was between the 
BiO and SrO planes. Remarkably enough, the DFT 
shows the presence of a narrow band of states near 
-1 eV induced by the dopant oxygen, as observed 
in experiment [1]. These states have 0(dopant)2pz 
character, possibly explaining their strong coupling 
to the STM tip states, but mix with Bi and 0(Bi) 
states as well. With the results of these calcula- 
tions, one may now attempt to correlate changes in 
gap observed by STS with changes in fundamental 
electronic structure parameters like t, t', electron- 
phonon coupling constants and Madelung energies. 

In the remainder of the paper we provide further 
evidence in favor of the hypothesis that the oxygen 
dopant atoms modulate the pair interaction locally 
rather than to act solely as potential scatterers. For 
this purpose we calculate the correlation between 
the local density of state (LDOS) and the modu- 
lated pair interaction/disorder potential and com- 
pare it with experimental results for the correlation 
of the LDOS with the positions of the oxygen dopant 
atoms as found by McElroy et al. [1] . 

2. Models of inhomogeneity 

Our calculations are based on the standard mean 
field Hamiltonian for a singlet d-wave superconduc- 
tor 

i]<y (ij) 

where we keep only nearest t and next-nearest t' = 
— 0.3t neighbor hopping, with /i = — t to model the 
Fermi surface of BSCCO near optimal doping. In the 
following we consider two models of how the oxy- 
gen dopant atoms create inhomogeneity in the lo- 
cal electronic structure. In model (i) we follow the 
conventional assumption that each oxygen dopant 
acts as a potential scatterer. The sum of all dopants 
gives rise to a screened Coulomb potential Vi = Vb/i 
at site i where fi = exp(-^is/'^)/'ris, with r^s 
the distance from dopant s to the lattice site i in 
the plane, in units of ^/2a, where a is the Cu-Cu 
distance. The pair interaction which enters the 
self-consistency condition for the nearest neighbor 
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d-wave order parameter Ay = 9ij{^i'\Cji — CjiCii) 
is assumed constant, i.e. gij — go in model (i). In 
model (ii) we assume that each oxygen dopant lo- 
cally modulates the pair interaction, which we model 
Qij = go + Sg{fi + fj)/2 with fi as defined above. 
In model (ii) the oxygen dopants exert no conven- 
tional Coulomb potential, i.e., Vi = 0. Since poten- 
tial scattering enters in the T3-channel in Nambu 
notation we will also call model (i) the T3-scattering 
model. Analogously we will term model (ii) the ti- 
scattering model since order parameter variations 
enter in the ri-channel in Nambu notation. 



3. Correlation between LDOS and dopants 

In Fig. 1 the LDOS is shown for both the ri- and 
the r3-scattering model along typical line cuts. The 
LDOS has been obtained by solving self-consistently 
the Bogoliubov-de Gennes equations resulting from 
Eq. (1) on a 80 X 80 lattice rotated by 45 degrees 
with respect to the Cu-Cu bonds, i.e., for a system 
containing 2 x 80 x 80 sites. Panel (a) shows the 
LDOS for the ri -scattering model. In good agree- 
ment with experimental STS-data [12], the LDOS 
along this linecut shows a strong anticorrclation be- 
tween height of coherence peaks and gap magnitude. 
Furthermore, the LDOS is quite particle-hole sym- 
metric and fairly homogeneous at small energies, 
which is also in favorable agreement with experi- 
mental findings. Panels (b) and (c) show the LDOS 
for the Ts-scattering model in the limit of a smooth 
potential variations, which arise if the dopants are 
poorly screened, and in the limit of a "spiky" poten- 
tial, which arises when each dopant is well screened 
and exerts only a short ranged potential. Obviously 
the T3-model fails to reproduce the experimental ob- 
servations in both limits. In the smooth limit (b) 
the height of the coherence peaks is positively cor- 
related with the gap magnitude in striking contrast 
to the experimental observation [1]. In the "spiky" 
limit (c) gap modulations appear to be quite small, 
the LDOS lacks particle-hole symmetry and the ho- 
mogeneity at small energies is destroyed by the for- 
mation of sub-gap resonances. 

A more direct test for the proposed scattering 
models is to consider the spatial correlation func- 
tions Cpg {Cpv) between the LDOS Pi{oj) and the 
modulated pair interaction gi (impurity potential 
Vi) which track the O dopant positions in our model. 
For example, the spatial correlation function be- 
tween the LDOS Pi{uj) and the impurity potential 
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Fig. 1. LDOS along typical line cuts for (a) ri-model with 
Vz = 0.57, A = 0.5, Sg = t, (c) smooth ra-model with Tz = 2, 
A = 2, Vb = l-5t and (c) "spiky" T3-model with Tz = 0.57, 
A = 0.5, Vo = t. 

Vi is defined as 

r i , {p^{u:)-{p{u:))){V-{V)) 

where Ay = ~ (^))^ {V) is the average 

oiV. 

The correlation functions for the ti- and r3- 
models are shown in Fig. 2 in comparison with the 
experimental result (panel (a)). The experimental 
correlation function between the LDOS and the po- 
sition of the oxygen dopant atoms (a) is remarkably 
particlc-hole symmetric with maximal negative val- 
ues near ±35 meV and maximal positive values near 
±70 meV. This behavior of the correlation function 
reflects the enhancement of the gap close to the oxy- 
gen dopant atoms, i.e., the coherence peaks which 
correspond to additional spectral weight move to 
higher energies. At small energies the correlation 
almost disappears, which is a signature of the ho- 
mogeneity of the LDOS at small frequencies or of 
the presence of other impurities uncorrelated with 
the oxygen dopants (see below) . 

The dashed line in Fig. 2(b) shows the correla- 
tion function Cpg between the LDOS and the pair 
interaction gi = ga = go + Sgfi for the Ti-scattering 
model. In good agreement with the experimental 
curve (a) , the correlation function is quite particle- 
hole symmetric, and is negative for small energies 
and positive for large energies. In contradiction to 
experiment, however, the correlation function re- 
mains negative and almost constant at small ener- 
gies, in particular it does not vanish for energies 
close to zero. This seems to indicate that the Ldos 
tracks the local gap magnitude at small energies as 
one would expect for a smoothly varying pair poten- 
tial, i.e., the slope of the LDOS at small energies de- 
creases because the gap increases in regions of large 
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Fig. 2. Correlation function between LDOS and (a) position 
of the oxygen dopant atoms (experimental) (b) pair inter- 
action Qi = go + Sgfi (ri-model with = 0.57, A = 0.5, 
5g = t), (c) impurity potential Vi (smooth ra-model with 
= 2, A = 2, Vb = 1.5t) and (d) impurity potential Vi 
("spiky" rs-model with = 0.57, A = 0.5, Vq = t). 



pair interaction (close to the oxygen dopants). We 
believe that this discrepancy with the experimen- 
tal data is due to the fact that the real BSCCO- 
2212 system contains further sources of disorder be- 
sides the oxygen dopant atoms. In particular, about 
0.2% in-plane native defects are typically imaged on 
a BSCCO-2212 surface. Due to the resonances they 
generate near zero bias [13] they are generally con- 
sidered to be strong pointlike scatterers. In order 
to mimic the real BSCCO-2212 system these addi- 
tional native defects have to be taken into account. 

The solid line in Fig. 2(b) shows the correlation 
function Cpg between the LDOS and the pair inter- 
action for a system where 0.25% strong pointlikc po- 
tential scatterers {V ~ hi) have been included in ad- 
dition to the dopants which only modulate the pair 
interaction. Obviously, the presence of a small num- 
ber of unitary scatterers affects the LDOS primarily 
at very small and at large energies. In particular it 
destroys the correlation near zero bias. The quali- 
tative agreement between the solid line in panel (b) 
and the experimental curve in panel (a) is almost 
perfect. Only the quantitative agreement is not as 
satisfactory because the correlations in our model 
calculation are by a factor of 2-3 larger than exper- 
imental observations. We speculate that this is due 
to additional disorder in the real BSCCO-2212 sys- 
tem which could e.g. stem from the random substi- 
tution of Bi at the Sr sites [2] . 

For comparison the correlation functions Cpv be- 



tween the LDOS and the impurity potential Vi for 
the T3-scattcring model arc shown in Fig. 2(c) and 
(d) in the smooth and "spiky" limit respectively. 
Both cases are in striking contrast to the experimen- 
tal data displayed in panel (a) . The most drastic dis- 
crepancy is the asymmetry of the correlation func- 
tion in the ra-model with respect to zero bias which 
simply reflects the particle-hole asymmetry of con- 
ventional potential scattering and is in clear contra- 
diction to the experimental observation. Only scat- 
tering at order-parameter variations can account for 
the remarkable particle-hole symmetry observed ex- 
perimentally. 

4. Conclusions 

By comparing the correlation between the LDOS 
and the oxygen dopants with experimental results 
we have provided further evidence in favor of a 
model where the nanoscale inhomogeneity in the 
local electronic structure of BSCCO-2212 is caused 
by a dopant modulated pair interaction. Our re- 
sults demonstrate also that the presence of native 
defects, which are imaged as scattering resonances 
near zero bias on typical BSCCO-2212 surfaces, 
reduces the correlation between the LDOS and the 
dopant atoms at small energies. 
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